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Dynamical mean field theory calculations are used to show that for stoichiometric, late transition- 
metal-oxides a critical variable for the Mott/charge- transfer transition is the number of d-electrons, 
which is determined by charge transfer from oxygen ions. Insulating behavior is found only for a very 
narrow range of d-occupancy, irrespective of the Coulomb repulsion. For representative Ni and Cu- 
based transition metal oxides, density functional band calculations yield d occupancies very far from 
the values needed to sustain an insulating phase, indicating either that density functional theory 
fails to predict the physical d occupancy or that the conventional Mott /charge-transfer concept does 
not solely account for the observed insulating behavior in these systems. 



I. INTRODUCTION 

The strong electronic correlations characteristic of 
transit ion metal oxides pose a problem of long-standing 
interest!^ and great current importance^'^. In transi- 
tion metal oxides the important low energy electronic 
states are derived from the transition metal d-orbitals. 
In a simple ionic picture the ground state corresponds 
to a d" configuration for each transition metal ion, with 
d-electron number n depending on material. The correla- 
tion energy cost involved in changing the d-occupancy is 
conventionally denoted by U, is defined a,sU — E[d^^^] + 
i?[c?""^] — 2E[d"'] and is typically large. In a ground- 
breaking paper MotiP^ argued if U were large enough 
relative to the bandwidth then electron repulsion would 
prohibit conduction, at least at integer valence. 

Mott's concept was refined in an important way by 
Zaanen, Sawatzky and AllerP who noted that correlated 
electron materials of interest typically contain ligands 
(for example the O atoms in transition metal oxides) 
which are close enough in energy to the chemical poten- 
tial that one must also consider the transfer of an electron 
from a ligand to a transition metal site. Ref. |5| defined 
the charge transfer energy Ect = E[d^^^] — E[d^] — Sp 
with —Ep the energy needed to create a hole in a ligand or- 
bital and demonstrated that if U were large but Ect < U 
then the physics would be controlled by Ect with U play- 
ing a relatively minor role. Implicit in this argument is 
a neglect of correlations on the oxygen sites; the justifi- 
cation is that the number of ligands is large enough and 
the hole density on the ligands small enough that con- 
figurations involving doubly occupied oxygen sites may 
be neglected. The resulting 'charge transfer insulator' 
paradigm is central to the conventional understanding 
of the physics of transition metal oxides. 

While density functional band theory (DFTjP is the 



work-horse of materials science, theoretical study of 
charge transfer insulators and other correlated electron 
materials requires methods which include additional cor- 
relations. Such "beyond DFT" methodologies require 
the identification of a set of orbitals whose correlations 
are to be treated more accurately and a prescription 
for embedding the correlated orbitals into the full elec- 
tronic structure. Two widely used examples are the 
"DFT-hU"G2l and "DFT-HDMFT"!2l31 approaches. To 
date, when these methods are applied to transition metal 
oxides the correlated orbitals are chosen to be the tran- 
sition metal d-levels. An important aspect of the the 
DFT+U and DFT-hDMFT approaches is the 'double 
counting correction t^^^, which removes correlation con- 
tributions from the DFT single-particle energies. The 
double-counting correction plays a particularly crucial 
role in charge transfer insulators because it affects the 
value of the charge transfer energy which is crucial to 
the physics. However, determining the double counting 
correction has proven problematic. Different prescrip- 
tions have been proposecP^t^^, but there is no consensus 
on what is the correct procedure. 

In this paper we argue that the double counting prob- 
lem should be recast as a problem of determining the 
occupancy of the correlated orbitals. Expressing the 
problem in this way reveals a fundamental and appar- 
ently previously unnoticed dichotomy in the materials 
theory of late transition metal oxides such as the high Tc 
cuprates: either the d-occupancy is very different from 
that predicted by density functional band theory or the 
materials are very far from the Mott-Zaanen-Sawatzky- 
AUen metal-charge transfer insulator phase boundary. 
The first alternative would be surprising because the 
present implementations of density functional band the- 
ory, while not exact, are generally believed to provide 
reasonable good approximations to the electron density. 
The second alternative challenges the currently prevail- 
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ing theoretical paradigm that the physics of cuprates is 
controhed by the Mott/charge-transfer transition. Our 
results also highlight a difficulty with the "DFT+DMFT" 
approach to transition metal oxides'^"*: the physics pre- 
dicted by this method is to a large degree controlled by 
the choice made for the double-counting correction. 



II. THEORETICAL APPROACH 

The approach we tak e is b ased on the "p-d" model, 
accepted in the literatur^^HESl a reliable description of 
transition metal oxides. The model includes correlated 
("d") orbitals with appropriate on-site interactions, lig- 
and (oxygen) orbitals (assumed uncorrelated because the 
hole density is small) and the relevant hybridizations. It 
may be written as 

H = Hd+ Hkyb + Hiigand (1) 

The Hamiltonian for the correlated ("d") subspace is 

Hd^Y. ^d4aAa. + (tV^ -l)+Hj (2) 

with d\^^ creating an electron in orbital a and spin a 

on site z, = X^ao- '^aa'^iacr the uumbcr operator for 
d electrons on site i and Hj representing the additional 
'Hunds-rule' interactions which give the multiplet struc- 
ture at fixed Nd (for explicit forms see e.g. Ref. f2l). The 
localized nature of the d orbitals means that Hj interac- 
tion terms are given by measured gas-phase values, which 
agree with calculationJSl, xhe 'charging energy' U is 
renormalized significantly by screening and is less well es- 
tablished theoretically although interesting results have 
appearecJ^. In this paper we explore a range of U but we 
emphasize that in the charge transfer regime of principal 
interest here the precise value of U is not important as 
long as it is large. The model neglects interactions on the 
oxygen sites, on the grounds that the number of holes per 
O site is low. While the question of the importance of the 
oxygen U should be reexamined with modern methods, 
this issue is beyond the scope of this paper. 

Hiigand and Hfiyi, describe the embedding of the cor- 
related orbitals into the broader electronic structure and 
it is generally accepted'^ ^ that reasonable estimates can 
be obtained from density functional band calculations, 
for example via the use of maximally localized W annier 
function techniques with a suitable energy window.^^^ 

The crucial model parameter is the difference A = 
Ep — Ed between an average ligand on-site energy Sp and 
the orbitally averaged d level energy Ed = Yla^dl^orb- 
Solution of the many-body model yields a physical charge 
transfer energy Ect which depends on A and on the 
many-body interactions. A direct fit of the quadratic 
parts of Eq. [T] to a density functional calculation yields 
Eqt = A; this value includes the density functional 
approximation to the interaction contributions, which 
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FIG. 1: (Color online) Metal-insulator phase diagram of 
model of high Tc cuprates (left panel) and rare-earth nick- 
elates (right panel) in plane of d — d on-site Coulomb inter- 
action strength U and d — p charge-transfer energy A. Filled 
symbols: limit Ac2 of stability of metallic (M) phase. Open 
symbols: limit Ad of stability of insulating (I) phase 



are also incorporated in the solution of the many-body 
model. Thus it has been recognized for many years that 
to obtain the Ed in Eq. [T] a "double counting correction" 
must be applied to the band theory estimatd^. Differ- 
ent prescriptions for the double counting correction have 
been proposecJlSlQll^ but no clear consensus has emerged. 
We investigate a range of A here. 

To solve Eq. [T] we use the single-site dynamical mean 
field approximatiorfiS! with the 'matrix' hybridization- 
expansion solver introduced in 1171 Some results were 
cross-checked by the exact diagonalization methocP^. 
The single-site approximation is generally accepted as 
containing the essence of the Mott and charge-transfer 
metal-insulator transition physics and is widely used in 
"DFT-^DMFT' studies of real materialfpHI. Muffi-site 
("cluster") methods, which are prohibitively computa- 
tionally expensive for multiband models with rotationally 
invariant interactions, provide a picture which is different 
in detaiP^ but confirm that the single-site approximation 
captures the important high energy physics of the metal- 
insulator transition. In particular a reasonable working 
definition of whether or not a material is "strongly corre- 
lated" is whether the stoichiometric compound is on the 
metallic or ins ulating side of the single-site DMFT phase 
boundarjIiSIini. 

We compute the behavior of two representative classes 
of charge-transfer materials, namely the layered high-Tc 
cuprates such as La2C'uOi and the pseudocubic rare- 
earth nickelates such as LaNiO^. For La2CuOi we use 
the standard "three-band" model with one d orbital per 
unit cell (representing the Cu d^^-y-i orbital) and two 
p orbitals (representing the orbitals on the two pla- 
nar oxygen ions). The Cu — O and O ~ O hoppings that 
make up H^v^+Hugand are obtained from band structure 
calculationJSIIl], There is general agreement on the Cud- 
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Op^ hopping parameter tpd ~ 1.6eV; different oxygen- 
oxygen lioppings are obtained by different workers but 
we have found that the behavior of the model does not 
in fact depend much on which value is choserP^. Because 
we consider only one d orbital in the cuprate case the 
only interaction is the Hubbard U. For LaNiOs we used 
a model with two d orbitals per unit cell, representing 
the — and — states, as well as the Op^ or- 
bitals on all three nearest-neighbor oxygen ions. We took 
the hybridization and ligand parameters from the GGA 
band calculations of Ref [53] (which we have verified are 
consistent with those arising from fitting the DFT results 
with maximally localized Wannier functions defined over 
the energy region 361^ above to 7eV below the chemical 
potential containing the p and d bands) . We included the 
full Slater-Kanamori interactions in the Cg manifolcP and 
set the exchange term J to 0.5 eV. In each case we fix 
the total number of electrons to correspond to the stoi- 
chiometric compound: one hole per unit cell in cuprates 
and three holes per unit cell in the nickclatcs. 



III. RESULTS 
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FIG. 2; (Color online) Metal-insulator phase diagram indicat- 
ing limit of stability of metallic phase Nd^c2 (solid symbols), 
limit of stability of insulating phase Nd,ci (open symbols) 
along with vertical lines indicating band theory estimates of 
Nd for model of high Tc cuprates (black and red squares, and 
black dashed line) and nickelates (blue and green triangles, 
and blue dotted line). 



The two panels of Fig. [l] plot the metal insulator 
phase diagrams calculated for the two models in terms of 
the traditional variables of interaction strength U and d- 
p energy difference A. The left panel agrees with r esults 
previously presented in the literature for cuprate^^aia^j 
The corresponding results for nickelates shown in the 
right panel have not been previously presented. In single- 
site dynamical mean field theory the metal insulator tran- 
sition is first order and is characterized by two spinodal 
lines: Uci{Aci) where the insulating phase loses stabil- 
ity and a slightly different value C/c2(Ac2) at which the 
metallic phase loses stabilitjH^. We present both lines, 
however we note that the Uc2 phase boundary has a 
strong temperature dependence which is computationally 
very expensive to capture. 

In Fig. [2] we replot the results in the plane of U and 
d-occupancy Nd with A as an implicit parameter. In the 
cuprates all of the d orbitals are filled except c?^2_y2 and 
we measure N^ with respect to the d^ core s < Nj, < 2 
with the physical range being 1 < Nd < 2. In the 
nickelates all d-orbitals are filled except the d^-i^y-i and 
dsz^-r^ and we measure Nd with respect to the d^ core 
so < Nd < 4 with the physical range being 1 < Nd < 4. 
The replotted results are seen to have a remarkably sim- 
ple structure. At small U the phase boundaries extrap- 
olate linearly to [/ = 0, Nd = 1 because as Nd — 1 the 
d bands become far removed from the p bands so both 
the d bandwidth and Nd — 1 ^ 1/A. However, for physi- 
cally relevant U values the phase boundary is essentially 
vertical, and occurs at Nd ^ 1.15 — 1.25 in the two di- 
mensional (cuprate) case and at Nd ^ 1.25 — 1.35 in the 
three dimensional (nickelate) case. 

The simplicity of the U —Nd phase diagram, along with 
the fact that the bare charge transfer energy A is a purely 



theoretically defined quantity which cannot be directly 
measured, identifies Nd as the critical variable for plac- 
ing transition metal oxides on the metal-insulator phase 
diagram. Further evidence that Nd is the important pa- 
rameter can be found in the electron spectral function 
(many-body density of states), obtained by maximum 
entropy analytical continuation from our imaginary time 
QMC data using the approach in Ref. 155 and shown for 
the nickelate case in Fig. [3] for several different combina- 
tions of interaction strength and charge transfer energy. 
The similarity of spectra with similar Nd but different A 
and U is evident. 



IV. CONCLUSIONS 

Nd has a precise meaning in the model defined by Eq. [T] 
but different definitions are available in band theory, in- 
cluding computing the charge density from orbitally pro- 
jected wave functions within a sphere of a given radius, 
and constructing maximally localized Wannier functions 
within some appropriately chosen energy window^^^. 
We have used literature values for the cuprate^^^'^il and 
obtained Nd for LaNiO^ from the maximally localized 
Wannier function construction defined above; the results 
agree within 10% with the Nd obtained from the sphere 
construction and reported in Ref. 1231 The results are 
Nd ~ 1.65 for cuprates (indicated by dashed line in 
Fig. [2]) and Nd ^ 2 for nickelates (indicated by dotted 
line in Figs. [2]). For both materials, the Nd values we ob- 
tain are so much higher than the critical values needed 
to drive the metal-insulator transition in single-site dy- 
namical mean field theory that any ambiguities in the 
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calculational definition of Nd are not relevant. 

However, the density functional results may not be 
definitive. While the approximate DFT implementations 
in use today are believed to produce accurate estimates 
of many electronic propertied, there is no solid-state 
benchmark for the quantitative reliability of the results 
of present state of the art functionals and some experi- 
mental results indicate that density functional methods 
overestimate covalencj^^fil. Thus it might be that density 
functional results for are substantially in error. This 
view is implicitly adopted in some of the DFT-I-DMFT 
literature. For example Ref.^ models La2Cu04 using a 
double counting correction which shifts from ~ 1.65 
to - 1.1. 

Experimental determination of the d valence is there- 
fore desirable. Indirect measurements, for example of the 
transferred hyperfine coupling^^ZI qj- the magnetic form 
factor in ordered states'^ have indicated very strong co- 
valency effects in cuprates. The estimates given in these 
papers are Nj, ^ 1.6, comparable to the band theory 
values. A more direct approach would be to detect the 
density of holes on oxygen via resonant X-ray scattering 
techniques as in Refs. I7|29l As indicated in Fig. [Sj spec- 
troscopic measurements of the relative position of p and 
d features in the many-body density of states would also 
be revealing. One particularly informative experiment 
would be a measurement of changes in the density across 
the ReNiOj, series, to reveal the change in O-hole den- 
sity. A systematic evaluation of the reliability of DFT 
estimates of transition metal-ligand covalency is also in- 
dicated. 

However, if the DFT results for Nd are reasonable, then 
associating the observed behavior of transition metal ox- 
ides with the Mott/charge- transfer insulator paradigm 
becomes problematic. La2CuOi is observed to be an in- 
sulator in experiment, and therefore if the d-occupancy 
predicted by DFT is not seriously in error, physics be- 
yond simple Mott/charge-transfer picture must be in- 
voked to drive the insulating state. For example, at mod- 
erate correlations a crucial effect of intermediate-ranged 
antiferromagnetic correlations has been found in clus- 
ter dynamical mean field and other approximations .f^SlEll 
LaNiOs is a metal in experiment, consistent with our 
GGA-hDMFT calculations based on the GGA Nd- How- 
ever, replacing La by other rare earths drives a metal- 
insulator transition which is sometimes associated with 
the Mott transition. Changing the rare earth distorts the 
lattice, and it is possible that this distortion might change 
Nd, but we have verified that within GGA band theory 
the d-occupancy is nearly the same for LuNiO^ in the ex- 
perimental structure as it is for LaNiO^. The results of 
this paper suggest that the observed insulating behavior 
is not solely determined by the Mott/charge-transfer phe- 
nomenon, but is instead associated with other physics, for 
example with the observed disproportionation in Ni — O 
bond lengths. GGA -|-U calculations (see Ref. [34] and to 
be published) confirm this association. 

To summarize, in this paper we have identified the 
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FIG. 3: (Color online) Local spectral functions for the nick- 
elates. Solid blue lines show the d-density of states and the 
dashed red lines show the p-density of states. The parameters 
for different panels are indicated. 



d-occupancy as a key variable for placing transition 
metal oxides in the Mott-Zaanen-Sawatzky- Allen metal- 
insulator phase diagram. We show that the conven- 
tional metal-insulator transition only occurs over a nar- 
row range of d-occupancies, that a transition will not oc- 
cur for sufficiently large d-occupancies, and that standard 
band theory techiques predict d occupancies which are 
deep in the metallic regime for La2CuOii and LaNiO^. 
While the notion of 'orbital' is not precisely definable in 
the solid-state context, ambiguities in the definition lead 
to uncertainties which are small on the scale of the differ- 
ence between the DFT Nd and the phase boundary. Thus 
the observed insulating behavior of La2CuOi indicates 
either that density functional band theory methods pro- 
vide remarkably erroneous estimates of Nd or that some 
physics beyond the Mott /charge-transfer paradigm is re- 
quired to explain the insulating behavior of this material. 
Our results have implications for 'beyond-DFT' theories 
of correlated electron materials. The beyond-DFT the- 
ories available to date typically require a definition of 
correlated orbitals which implies precisely defined orbital 
occupancies, and must introduce double-counting correc- 
tions. As we have shown, it is these double-counting cor- 
rections which fix the d-occupancy and therefore dictate 
the resulting physical behavior. 
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